Discrete Fourier Transform in Nanostructures using Scattering 
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In this paper we show that the discrete Fourier transform can be performed by scattering a coherent 
particle or laser beam off a two-dimensional potential that has the shape of rings or peaks. After 
encoding the initial vector into the two-dimensional potential, the Fourier-transformed vector can 
be read out by detectors surrounding the potential. The wavelength of the laser beam determines 
the necessary accuracy of the 2D potential, which makes our method very fault-tolerant. 
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I. INTRODUCTION 

The fast developing fields of quantum computing and 
quantum information processing have attracted much in- 
terest recently 1-6 . Most of the papers in the field of quan- 
tum computing, such as Shor's and Grover's about fac- 
toring numbers 7 and searching a database 8 , use many- 
particle physics (i.e. qubits) for their implementation. 
However, it is also possible to implement quantum in- 
formation processes in the so-called unary representation 
of dimension TV, which consists of a statistical ensemble 
of non-interacting TV-level systems that can be described 
by a single-particle Hamiltonian. So only interference, 
but no entanglement is needed for the implementation of 
quantum algorithms in the unary representation. It was 
already shown that the Grover algorithm 9 can be per- 
formed in a single shot using the unary representation. 
While in Ref. 10 the quantum database search has been 
done experimentally with atoms in a beam that serve as 
non-interacting AT- level systems, Refs. 11 and 12 show 
theoretically that the large spin of molecular magnets 
and nuclear spins in semiconductors can be used as non- 
interacting TV-level systems, respectively. Here we pro- 
pose a method to perform the classical discrete Fourier 
transform (DFT) in the time ensemble of quantum or 
classical scattering events, where the scattering poten- 
tial provides the classical TV-level system for the read-in 
of information and the ./V detectors provide the classical 
TV-level system for the read-out of the information. In 
the case of quantum scattering a single-particle Hamil- 
tonian describes the information processing. In contrast 
to Shor's algorithm our approach requires an exponential 
overhead of the hardware. However, once our device is 
built the DFT can be computed very efficiently, i.e. in 
a single shot. In our method we propose to scatter a 
coherent matter or laser beam off a two-dimensional po- 
tential consisting of TV concentric rings or TV peaks, the 
amplitudes of which represent the elements of a vector in 
an TV-dimensional vector space where the information is 
encoded. The Fourier transformation of this vector can 
then be read out immediately by detectors surrounding 



the 2D potential. Our method is closely related to the 
general field of Fourier optics 13-15 , where a diffraction 
pattern is Fourier transformed into a desired image pat- 
tern by means of laser scattering. We go a step further 
in the sense that we describe how to use a 2D scattering 
potential to perform the DFT efficiently. The concentric 
rings and peaks can be produced by charge or spin distri- 
butions. The calculation time t c is roughly given by the 
diameter of the detector ring (e.g. 1 mm) divided by the 



light velocity c = 3 x 10 m/s, i.e. t c 



10" 



s. Using 



for example a laser with a wavelength of 500 nm and a 
spot radius of 100 /nm enables the Fourier transformation 
of a function with 100 sampling points. The larger the 
spot size and the smaller the wavelength of the laser, the 
more sampling points can be achieved. The advantage 
of this device would be that it can be used already for 
example as DFT coprocessor to speed up present com- 
puters operating at room temperature. Rough estimates 
show a possible speedup of 10 6 or more (see below). 

In Sec. II we give an overview on the DFT and the Fast 
Fourier Transform (FFT). Sec. Ill and Sec. IV describe 
the quantum-mechanical coherent matter scattering and 
the classical laser beam scattering off a specially designed 
2D potential that is used for the implementation of the 
DFT. 



II. DISCRETE FOURIER TRANSFORMATION 

The discrete Fourier transformation (DFT) has a wide 
range of applications in the field of digital signal process- 
ing, such as spectral analysis and filtering. The DFT 
enables the possibility to extract the period of a func- 
tion, which can be used to factor numbers. Thus finding 
efficient ways to compute the DFT is the key to the ef- 
ficient factorization of numbers. The calculation time of 
a DFT can be greatly reduced by the classical algorithm 
called Fast Fourier Transform (FFT), which gets rid of 
the redundancies found in the DFT. In a similar way 
Shor's quantum algorithm reduces the calculation time 
of the discrete quantum Fourier transformation (DQFT) 
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by means of an efficient quantum circuit 7 that uses qubits 
as information carriers. In our proposal the superposi- 
tion principle of quantum or classical wavefunctions is 
used to perform an TV-dimensional DFT with N shots 
or even a single shot in one scattering device. The FFT 
provides then the possiblity to link several scattering de- 
vices in parallel. This means that the speedup factor N 
can be multiplied by the number of the parallel scattering 
devices. 

First it is instructive to give the definition of a DFT. 
When we calculate the DFT of a function f(z) that maps 
complex numbers on complex numbers, we can first dis- 
cretize the function f(z) in M total sampling points 
z 1 ,z 2 ,...,z M , leading to M values /i = f(z 1 ),f 2 = 
f(z 2 ),...,f M = f(z M )- Then the DFT of / is given 
by 

M-l 

F k =J2 e 2 ^ klM ! r (1) 
j=o 

The matrix elements of a DFT consist of the irreducible 
representation of the cyclic group produced by a rotation 
R-2tt/m about the angle 2n/M, where M is the dimension 
of the vector space, which is equal to the number of sam- 
pling points. For example a DFT of dimension M — 4 is 
represented by 
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It will turn out that the exponential overhead of the hard- 
ware imposes a physical limit to the number of sampling 
points N that can be processed by our DFT device. How- 
ever, the Danielson-Lanczos Lemma 16 , which is widely 
used for the FFT algorithm, states that the DFT for M 
sampling points can be separated into the sum over the 
DFT for the M/2 even sampling points and the DFT for 
the M/2 odd sampling points, i.e. 

M/2-1 M/2-1 

F k = Y, e 2 * ik W M f 2j + Y e 2mk ^ +1 ^ M f 2j+1 . 

3=0 j=0 

(3) 

Usually, the FFT uses this Danielson-Lanczos Lemma re- 
cursively down to Fourier transforms of length 2°. There- 
fore it is essential that M = 2™, where m is a positive 
integer. As we want to use 2D scattering potentials with 
real amplitudes, we can separate the real and imaginary 
part of the DFT: 

M-l M-l 

F k =Y e 2 ^ k ' M g 3 + i Y e 2 ^ k ' M h^ (4) 

j=0 j=0 

where f(zj) — g(zj) + ih(zj) with real functions g(zj) 
and h(zj). 



Since our device can perform the DFT efficiently for 
N sampling points, we can stop the recursion at N = 2 n 
sampling points, where n is a positive integer that is 
smaller than m. Using our device as coprocessor of an 
existing computer, the remaining work of the main pro- 
cessor consists in combining the M/N partial Fourier 
transforms, which takes of order log 2 (M/N)M/N — 
(to — n)M/N operations. Thus, compared with initially 
necessary Mlog 2 (M) operations a speedup of at least a 
factor N is achieved. In this estimate of the speedup, 
the necessary time to calculate the DFT for the N sam- 
pling points is assumed to be of the order of 10 ps, i.e. 
negligible. Alternatively, one could say that an existing 
computer equipped with our device is capable of factor- 
ing numbers that are N times larger than the largest 
nowadays factorizable number. If M/N DFT devices are 
linked in parallel, a speedup of factor M is achieved. 




FIG. 1. 2D scattering 

potential V(x') = (1 + sin a) 53 =1 ajS(r' — r^/Vr'j which 
performs the DFT in N shots in the case of elastic scattering 
and in a single shot in the case of inelastic scattering. 

III. COHERENT MATTER SCATTERING 

Now we describe the various possible implementations 
of our DFT device, in which a coherent matter or laser 
beam is scattered off a 2D potential. This potential con- 
sists of radial 5- functions or <5-peaks on a straight line, the 
amplitudes of which both sample the points of the real 
or imaginary part of our function f(z) and represent the 
amplitudes of our basis vectors. Let v = Y] i a j w j be 
the input vector, where Wj are the basis vectors with am- 
plitudes a,j — g(zj) or dj — h(zj). Then our 2D potential 

is of the form U(x) oc Ylj=i a j$( x ~ x j)- The read-out is 
done by detectors surrounding this potential. While the 
coherent matter beam needs a high precision of the po- 
tential, the laser beam provides a high fault-tolerance of 
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the potential since faults that are smaller than the wave- 
length of the laser beam remain undetected. We give 
the derivation for both the coherent matter and the laser 
beam scattering. We first start with the coherent matter 
beam. We use the nomenclature of Ref. 17. First we cal- 
culate the 2D scattering amplitude, from which we can 
infer the precise shape of the 2D potential required for 
the implementation of the DFT. For this we have to start 
from the well-known Lippman-Schwinger equation 17 



x 



(±) 



(x|0) + |j J dVG±(x,x')(x'|T/|v (±) 



(5) 



where \(j>) is the incoming wave, V is the 2D scattering 
potential, and 

is the Green's function, which also solves the Hclmholtz 
equation (V 2 + fc 2 )G±(x,x') = 6 (2 ^(x - x'). For the 
evaluation of the Green's function we transform first to 
the momentum representation, i.e. 



G ± (x,xO = - 



d 2 p' 



4p'-( x - x ' 



(2nh) 2 E 



2m 31 16 



(7) 



where we have used the 2D overlap (x|p) = e lp ' x ^ h /2nh. 
Inserting E = h 2 k 2 /2m 7 p' = hQ, and going to polar 
coordinates (Q,<j>), we get 



G±(x,x 



(27T) 



QdQ 



QdQ 



giQ|x— x'| cos0 

' k 2 -Q 2 ±ie 



2 3/2 



cos(Q|x-x'|-f) 
k 2 -Q 2 ± it 
I e ±«( fc | x - x '|-f) 



4^ y/2k |x - x' 



(8) 



For the approximation we have used the integral rep- 
resentation of the Bessel function Jo(Q|x — x'|) = 

h t e iQ l x ~ x 'l « vS-^ C0S(g |X X ' h ^ 
which holds in the asymptotic regime Q |x — x'| >> 1, 
i.e. the local 2D potential (x'| V |x") = F(x')<5 (2) (x'-x") 
and the detectors at x are far apart from each other. 
Then we obtain with Eq. (5) 
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where r := Ixl, r' := |x' 
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k' := fcA, and 



/(k',k) 



2tt 2m 
4~ J 2 



yrfVe- 4k '- x V(x')(x'|v- (+) 



(10) 



is the 2D scattering amplitude. The amplitude of the 
2D spherical wave decays as 1/ ' \fr. This is in agreement 
with the topological argument that the circumference of a 
circle is 2-kt and that the probability has to be conserved 
on the circle. In first-order Born approximation, where 

|/(k',k)| <C Vkr, we can insert (x' = e 2n into 

Eq. (10), and we define q := k' — k as the scattering 
wavevector. Eq. (10) becomes then 

/«(k',k) = -—j=—j d 2 x'e^ c °^F(x'), (11) 

where ip is the angle between q and x'. 




FIG. 2. If the scattering is elastic, the DFT can be per- 
formed with N shots. The 2D potential shown in Fig. 1 is 
located in the center of the ring. The length of the lines be- 
fore and after the scattering are proportional to the length of 
the wavevectors k and k', respectively. 

For the implementation of the DFT we start first with 
a 2D potential of the form V(x') = V(r')(l + cos^). If 
the scattering is elastic, we choose ip = n/2 — 9/2. If 
the material allows inelastic scattering, we choose ip = 
n/2 — 9. It will turn out that for this potential q has 
to point always into the same direction. For example q 
can be always perpendicular to the x-axis. Then we can 
define a as the angle between the x-axis and x', which 
turns our potential into V(x') = V(r')(l + sin a). This 
2D potential is shown in Fig. 1. Inserting this potential 
into Eq. (11) leads to 

/M (k ',k) = -^|2 J r'dr' [Jo(qr>)+iJi{qr')]V(r>), 

(12) 



which in the asymptotic limit r' 3> 1/q becomes 

/«(k',k) = -2^=|? / V?dr'e^ r '-Vv(r>). (13) 
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By choosing V(r') = YljLi a j^( r> ~ r j)/Vr* we can en- 
code the information into the N amplitudes a,j and per- 
form a DFT 

i 9 N 

' (I> < k ' k > = -^l ie - , \X> ,w <»» 

In order for it to be a discrete Fourier transformation 
(DFT), we have to discretize the variables. Let us 
make a first attempt: rj — j4?, j — 1, 2, . . . , N, and 
q v = v=f-, v = 1, 2, . . . , N, with v referring to the ^th 
detector. The smallest optically resolvable distance be- 
tween the radii rj is A/2. So it is safe to choose I = NX, 
where A = 2ir/k is the deBroglie wavelength of the co- 
herent beam. Then we obtain jt = 2 sin | for the elastic 
scattering and ■£> = tan for the inelastic scattering (see 
below) . The radii of the rings are multiples of the wave- 
length A, i.e. Tj = jX. There is still the requirement 
Tj ^> l/q v that has to be satisfied. Therefore we have to 
choose j = N + 1, N + 2, . . . , 2N, which still produces a 
DFT, since the DFT is 27r-periodic. 




FIG. 3. The 2D scattering potential shown in Fig. 1 is lo- 
cated at the center of the ring. The length of the lines be- 
fore and after the scattering are proportional to the length 
of the wavevectors k and k', respectively. The DFT can be 
performed with a single shot if the scattering off the 2D scat- 
tering potential is inelastic. 

Our 2D potential V(x.') has the 1 + cos<^ dependence, 
so that we obtain Jo and J\ Besselfunctions in Eq. (12). 
Since we do not want to adjust the orientation of V(x') 



according to the scattering vector q, we have to keep the 
direction of q constant in space. This requirement can be 
satisfied in two ways, which can be identified by elastic 
or inelastic scattering. In the case of elastic scattering k 
and k' have the same length. So q := |k — k'| = 2fcsin |, 
where 9 is the angle between k and k'. For every incom- 
ing beam of angle $i n = 9/2 the resulting vector element 
after the DFT can be read out in the detector at the out- 
going angle i? ou t = i?i n = 9/2. The number of incoming 
and outgoing angles is exactly N . Thus the DFT can be 
performed with TV shots, which is shown in Fig. 2. In the 
case of inelastic scattering, k and k' do not have the same 
length anymore. It is then possible to perform the DFT 
in a single shot. For this q := |k — k'| = fctan#. Since 
k' = n r k and k = k' cos 9, we obtain an angle-dependent 
refractive index n r = l/cos0, which is probably difficult 
to realize experimentally. At least in the case of a laser 
beam (see below) it may be possible to use photonic crys- 
tals that mimick n r = l/cos#. The inelastic scattering 
is shown in Fig. 3. 




FIG. 4. 2D scattering potential 

V(x') = 53^=i a j&i r ' ~ r j)S{ a ~ a j)/ r \ which performs the 
DFT in a single shot for elastic scattering. 



Instead of V(x') = (1 + sin a) Y^,f=i a j^( r ' ~ f^/Vr 7 
consisting of concentric rings, we can alternatively use 

peaks of the form V(x') = Ylj=i a j^( r ' ~ r j)3( a ~ a j)/ r ' \ 
which is shown in Fig. 4. After putting the peaks on a 
straight line {pt\ = oti = ■ ■ ■ = = 7r/2), we obtain 
from Eq. (11) 

/« (k', k) = --1= ^ Yl a ^ rj cos v " ' ( 15 ) 

v n j=l 

where a v = tp v + ir/2 is the angle between the scattering 
vector q„ and the x-axis. With this potential we can dis- 
cretize the variables similarly as before: Tj = j-h, j = 
1, 2, . . . , N, and q v cos ip v = u^f , v = 1, 2, . . . , N. The 
advantage of this method is that the peaks can be much 
easier produced experimentally than the rings. Also, we 
can start from the peak with j = 1, because we do not 
require Tj 3> l/q v anymore. There is also no scaling fac- 
tor 1/s/q^ needed for the detectors (see Eq. (14)) and 
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the scattering can be done elastically with a single shot 
(see Fig. 5) because the scattering vectors q„ do not need 
to point into the same direction. However, the potential 
made of peaks decreases as 1/r', and not as 1/Vr 7 as for 
the rings. 

For the experimental realization of the 2D scattering 
off a 2D potential one could use a 2D electron gas that is 
scattered off an electrostatic potential. Like in the case of 
the laser beam scattering (see next section), it would be 
also possible to use 3D scattering off a 2D potential, such 
as the scattering of spin-polarized neutrons off a spatial 
electronic spin distribution. 

IV. LASER BEAM SCATTERING 

Since coherent matter waves have a very short de- 
Broglie wavelength, the resolution of the 2D potential 
must be very high (< 1 A), which is beyond today's tech- 
nology. In order to overcome this difficulty, one can use 
3D scattering of a laser beam off a 2D potential, where 
the wavelength of the laser can be chosen according to 
the accuracy of the 2D potential. For this we have to 
start from Maxwell's equations. Following Ref. 18 we 
need to solve the Helmholtz equation 

(V 2 + fc 2 )D = -V x V x (D - e E) 

-!f^Vx (B- Mo H), (16) 

where D and B are the electrostatic and magnetic fields 
within the sample, E and H are the respective fields of 
the laser beam, e is the dielectric constant, fi is the 
permeability, and k — y/eonouj Jc is the wavevector of our 
laser beam. The result reads 18 

D = D<°) + A s e tkr /r, (17) 

where the scattering amplitude in Born approximation is 
given by 

Do 4W \ eo 

+ (nxe*)-(n xeo)^H, (18) 
Mo J 

where <5e(x') <C e and <5/i(x') <C Ho- Here n and e 
are the direction and the polarization of the incoming 
beam, n and e are the direction and the polarization 
of the outgoing beam, respectively. Let us restrict our- 
selves to variations in the dielectric constant Se(x'), i.e. 
S/j,(x') = 0. In analogy to the 2D potential for the coher- 
ent matter beam, we can use a dielectric variation of the 
ring form Se(x') = (1 + sina) X^jLi a j^( r> ~ r j)l^ ! or 01 

the peak form 5e(x') = a j^( r ' ~ r j)^( a ~ a j)/ r ' ■ 

Performing the integration in cylindrical coordinates, 
f d 3 x' = J r'dr'dipdz, for the peak form leads to 



e* • A (1) k 2 d 

D 47re t— f 

where d is the thickness of the sample. So 3D scattering 
off a 2D potential provides another means to perform the 
DFT. 

For the experimental implementation the spatial di- 
electric variation <5e(x') can be produced for example elec- 
trostatically by means of the Kerr effect. As another ex- 
ample one could use a crystal with spin-orbit interaction 
where an imbalance of spin-up and spin-down electrons 
leads to different dielectric constants for left and right 
circularly polarized light 22 . In this case a circularly po- 
larized laser beam is scattered off a spatial dielectric vari- 
ation <5e(x') that is created by a spatial electronic spin 
distribution. It is also worth to mention that the laser 
scattering can be done at room temperature. 

Now we estimate the amount of sampling points (=N) 
that can be used with today's technology. The state-of- 
the-art lasers of e.g. A = 500 nm wavelength produce an 
intensity of J=l GW/cm 2 at the sample. This implies a 
photon number density of n = <~ 10 10 photons per 
(100 /zm) 3 . If only 10~ 2 of the laser beam is scattered, 
about 10 4 photons remain at the detectors surrounding 
the sample at r = 1 mm. Within a radius of 100 /xm= 
2iVA we can insert N — 100 rings or peaks, i.e. use a 
basis of N = 100 sampling points. 




FIG. 5. The 2D scattering potential shown in Fig. 4 is lo- 
cated in the center of the ring. The length of the lines before 
and after the scattering are proportional to the length of the 
wavevectors k and k', respectively. The DFT can be per- 
formed with a single shot by using elastic scattering off the 
2D scattering potential. 

Until now we have taken a 2D variation of the dielec- 
tric constant <5e(x') = X^jLi a j3( r ' ~ r j)3( a ~ a j)/ r '- F° r 
experiments it is more realistic to assume a 2D varia- 
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tion of the dielectric constant <5e(x') = X^=i CLj^ w \r' — 
rj)5^ w \a — a.j)jwr' made of delta functions of width w, 
such as 5^ w \r' - rj) = w/n[{r' - rj) 2 + w 2 }. Then Eq. 
(19) becomes approximately 



e* • A 



(i) 



Dn 



■eo^ajj dr" 



k 2 d 
4ire 

k 2 de -q v \co Sl p v \w 
47T60 



iq v {?" +rj) cos <fi> 



+ W 2 



N 



i=i 



e y aje" 11 " 3 ^^, (20) 



where we have made the substitution r" = r' — rj. As 
long as w C l/2k = A/47T < l/q v | cos <p v \ , this expo- 
nential damping factor is negligible. Using a laser of 
wavelength A = 500 nm implies that w ^ 50 nm. It 
is also possible to scale up the proposed setup to longer 
wavelengths without changing the results. One could say 
that the required accuracy of the 2D potential is given 
by the wavelength of the laser. Therefore our scheme is 
very flexible and very fault-tolerant. 
Since jfc = 2 sin | for elastic scattering, the detectors are 
to be arranged within a scattering angle of 7r/3. The 
smallest angle between two detectors is m ; n = tt/3N — 
7r/300, which leads to a minimum distance of r9 m i n = 10 
/im between the detectors. 



V. CONCLUSION 

In conclusion, we have shown that by means of scat- 
tering a coherent matter or laser beam off a 2D potential 
it is possible to perform the DFT on TV sampling points. 
So the scattering is done quantum-mechanically or clas- 
sically. The 2D potential can be made of concentric rings 
or peaks on a straight line. The scattering off the peaks 
should be much easier to implement experimentally since 
it can be done elastically with a single shot. Note that 
the DFT is performed classically in both cases since we 
use a classical 2D potential to encode classical informa- 
tion. This provides the great advantage that our de- 
vice can already be used as DFT coprocessor for existing 
classical computers, without having to worry about the 
decoherence of quantum states, like in the case of quan- 
tum computing or spintronics devices 19-23 . In addition, 
M/N devices can be linked together to perform the DFT 
in parallel on M sampling points, thereby making use 
of the FFT. It should be possible with today's technol- 
ogy to build a DFT device with the capability to Fourier 
transform an TV = 1024-dimcnsional vector. If an exist- 
ing classical computer was equipped with M/N = 1024 
such DFT devices that are linked together, a speedup of 
a factor M = 1024/V - 10 6 could be achieved. 
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